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ABSTRACT 


Electronic Oscillators are studied as an application of the general 
method in the parameter plane of analyzing non-linear systems with self- 
excited oscillations. Computational methods are developed which avoid 
the lengthy calculations required by the general method for the second 
and third order systems which arise. The single parameter and multi- 
parameter sensitivity of the frequency characteristics are determined 
and it is demonstrated that these sensitivity measures may be used to de- 
termine the best configurations of networks to be used to minimize sensi- 
tivity to specific parameters. The method employed suggests application 
to the frequency characteristics of passive and active filters. 

A specific problem illustrating the failure of the © = O contour in 
the parameter plane to completely delineate regions of stable and un- 
stable solutions to the system differential equation is included in the 
appendix. No solution is given but suggestions as to how this apparent 


paradox arises, are made. 
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SENSITIVITY ANALYSIS OF SELF-EXCITED 
OSCILLATIONS IN NON-LINEAR SYSTEMS. 


AN APPLICATION TO ELECTRONIC OSCILLATORS. 


LE Introduction. 

The behavior of electronic oscillators has been described both 
qualitatively and quantitatively in many ways. However, as new models 
and methods in analysis are developed for the solution of new problems, 
it is worthwhile to test their application to older problems to see if 
a simplification or unification of theory exists. It is in this spirit 
that electronic oscillators are studied herein. They are considered as 
a class of non-linear active networks in which self-excited oscillations 
arise. 

The principle analytical tool to be used is the Parameter Plane 
Method. This application to electronic oscillators is, therefore, a 
practical application of general theory developed by Siljak[20] and Sil- 
jak and Stoic{21] for self-excited oscillations in non-linear systems. 

In practical oscillator design and application, the fundamental 
problem is amplitude and frequency stability with respect to changes in 
system parameters. Second to this may be the accurate control of fre- 
quency by deliberate variation of system parameters. This problem is 
analyzed by one parameter and multiparameter sensitivity analysis. In 
this respect, the sensitivity problem is approached in a historical se- 
quence. It is first considered as an afterthought to oscillator design 
and then possible design techniques are considered for use when the sensi- 
tivity is given as ana priori design specification. 

Several well known oscillator circuits of both RC and LC type are 


analyzed and sample problems are worked to show how the parameter plane 


representation and sensitivity analysis can be used to choose one top- 


ology over another to satisfy design specifications. 
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a. The Parameter Plane Method. 

Before the specific application of the parameter plane method is 
made, it is worthwhile to trace the development of the method as an 
algebraic technique in control system design and to examine the basic 
theory. 

Mitrovic{ 18] is generally credited with initiating the method in 
the early 1950's. His work was directed to finding the roots of a poly- 
nomial as functions of the coefficients of the first and zeroth order 
terms. (The role of finding the roots of a polynomial is a control 
engineering oriented problem, to say the least:) He found it possible 
to map contours of constant damping factor, €, and undamped natural fre- 
quency, w from the complex frequency domain into a plane with the co- 
efficients as ordinate and abcissa. However, considerable extension of 
theory was required in order to develop the parametric equations of 
suitable contours in terms of realistic system parameters which enter 
into many coefficients of a characteristic polynomial of a system. The 
method, as a graphical tool for finding roots and determining system 
Stability has much in common with the root locus technique, but with the 
facility to observe the movement of roots simultaneously varying with 
respect to two parameters rather than one. Attempts have been made to 
extend the concept into a parameter space to consider variation of roots 
with respect to more than two parameters, but interpretation of results 
rests on developing a greater facility with the parameter plane and the 
algebraic expressions associated with it. 

In the complex frequency domain (s-plane), a root position can be 
characterized by a set of cartesian or polar coordinates. Theory of 
the complex variable suggests the real and imaginary part of the com- 


plex variable, s, as Suitable cartesian coordinates and the practical 
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criteria of damping factor and undamped natural frequency borrowed from 
second order system analysis are suitable, if unconventional, polar co- 
ordinates. Figure 1 shows the relation of these coordinates. 

It is seen that 


Ss= "Ot jw 


or s= - w cos® + jw sin® 
n n 

but C= cos0 

hence sind = /1-¢2 


from which we have 
s=-o+ jw (1) 
se - wut t ju/i-@ (2) 
The general nth order characteristic equation of a system can be 
represented by 
= a,s* = 0 (3) 
ke O 
Mitrovic directed his effort to developing expressions for the paramet- 
ric equations of contours of constant % and w, in the coefficient plane, 
that is, the chosen parameters were the coefficients of the first and 
zeroth order terms in (3). The application to finding the roots of 
polynomials was direct, and in particular led to a single normalized 
graphical relation for all of the roots of any third degree polynomial. 
However, in application to control systems, it was seen that using co- 
efficients as parameters was less than satisfactory as a general method 
since the real system parameters appeared in combination in many co- 
efficients. 
Siljak[{19] generalized the method of Mitrovic to the case where the 
real parameters appeared in linear combination in any coefficient. If 


a and § are two parameters, in general, we might have 
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jw 


Figure l 


Coordinates of Point in s-Plane 
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ay = bo + cP + d (4) 


Subsequent generalization by Hollister[{23] and others has led to solu- 
tions for systems in which the parameters appear as products or to the 
second power. To date, this has been sufficient to allow the applica- 
tion of the method to a large number of problems in control system and 
network analysis. 

At the heart of the method is the simple fact that, since s is a 
complex variable, then any equation f(s) = O has a real and imaginary 
part which are independently zero. This gives rise to two simultaneous 
equations in any pair of parameters @ and § in terms of the coordinates 
of the s-plane, (¢ and On from (2) or o and w from (1)). 

Recent work by Siljak[20] has shown that this essentially linear 
technique can be extended to yield accurate results in the prediction of 
limit cycles and the stability of limit cycles in non-linear systems. 
The frequency of a limit cycle can be determined and, in some cases, the 
transient response of a non-linear system may be constructed from the 
parameter plane representation. The application is subject to limitae 
tions at the present time but is based on the conjecture that, if a lin- 
ear system can be represented by the collection of stationary roots of 
the characteristic polynomial, then a nonelinear system may be character- 
ized by a set of moving roots. Knowledge of the dynamics of root behav- 
ior should then yield knowledge of the system behavior in the time do- 
main. Of course, no definition exists of a "dynamic root" at this point 
in time, since the whole concept of a characteristic polynomial of a 
system is based on the essentially linear analysis technique of the La- 
place transform. Yet, the notion of moving roots is producing meaning-= 


ful results under certain restrictions. 
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The application of the parameter plane method to non-linear systems 
is just such an extension. Whereas, in a linear system we speak of a 
particular root system being characterized by an "M-point" which is fix- 
ed in the parameter plane by the choice of a parameter pair, we must now 
speak of a "dynamic M-point", the locus of the values of the parameters 
as dictated by the nature of the non-linearity. This is the technique 
to be used in the analysis of electronic oscillators. 

From Mitrovic's original idea has been developed a comprehensive 
analysis tool. From time to time one hears objections to the method as 
being unnecessarily complicated or even redundant in the face of other 
analytic methods. Granted, the method applied in general is best used 
where a general purpose computer is available to plot the parameter 
plane contours, yet, for many problems there are simplifications which 
show the utility of the method while avoiding the complexity. The appli- 
cation of the method to analysis of the sensitivity problem in frequency 


of oscillators is just such a case. 
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3. Electronic Oscillators. 

In oscillators producing a nearly sinusoidal wave-form, amplitude 
stabilization is achieved by the onset of some non-linear property of 
an element of the network. This could be the effect of saturation or 
the change in gain of the active element, and in some practical oscilla- 
tors it has been brought about by using a lamp as a non-linear resist- 
ance to affect the bias of a variable transconductance type vacuum tube 
to control the gain. In the usual analysis, it is assumed that the 
oscillation can be described in terms of the location of a pair of poles 
of the closed loop transfer function on the imaginary axis of the com- 
plex frequency plane. Such a delicate balance is highly unlikely to 
last very long in a truly linear system. The application of the para- 
meter plane analysis removes the necessity of the assumption, provides 
a new insight into the stability or instability of the oscillation 
which is properly a limit cycle, and then provides a validation for the 
above assumption for the case of oscillations with a nearly sinusoidal 


wave-form. 


Application of the Parameter Plane Method 
In this analysis, it is assumed that there is only one non-linear 

element. One of the two variable parameters is assigned to this non- 
linear property. If there are two non-linearities, a variable para- 
meter designation is assigned to each. If there are more than two such 
non-linearities, then assumptions must be made as to a representative 
single-value to be used for all but two non-linearities or extensions 
must be made into a "parameter space" where one would consider a "dy- 
namic M-surface" rather than the dynamic M-locus considered in the two- 


dimensional parameter plane. In the case of electronic oscillators, 
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there is usually only one significant non-linearity in the system, and 
this will be the case in all of the examples studied in this thesis. 
The second variable parameter is assigned to a system parameter of in- 
terest and all other parameters are fixed at some reasonable value. By 
mapping the stable and unstable regions of the parameter plane, and by 
superimposing the excursions of the parameter assigned to the non-lin- 
earity, it is possible to determine the existence of a limit cycle, to 
predict its stability and to predict the frequency. 

In the s-plane, the stable and unstable regions are bounded by the 
Ss = jw or (= 0 line. This contour can very easily be mapped into the 
parameter plane without recourse to the extensive general methods. By 
substituting s = jw into the characteristic equation of the oscillating 
aston’, the equation can be readily divided into a real and imaginary 
part, which equated separately to zero yield two simultaneous equations 
in the parameters, say a and §, as functions of the frequency w. These 
are the parametric equations of the contour ¢ = 0 or g = O in the para- 
meter plane. 

Since there is only one non-linearity in consideration, the M-locus 
is a straight line. The intersection of this line with the & = 0 con- 
tour indicates the existence of a limit cycle at the frequency given by 
the parameter w at the point of intersection. A qualitative argument 

hie is at this juncture that it is quite proper to object from a 
mathematical standpoint that this whole structure violates some very 
basic concepts. A characteristic equation arises from the Laplace Trans- 
formation of the system differential equation. If any coefficient or 
coefficients are time-varying (as they most assuredly will be in the case 
considered), they cannot be removed from the transformation integral and 
it is not correct to carry such coefficients into a "characteristic 
equation". I cannot answer such an objection, indeed anyone who can has 
a good handhold on the problem of analyzing non-linear systems as systems 
with moving roots. Yet this is the basis on which I will proceed and let 


the results themselves speak for the validity of the assumption and beg 
for a more complete mathematical justification! 
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Analysis of Stability of Limit Cycles 
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may then be given, based upon the nature of the non-linearity as to 
whether the limit cycle is stable or not. 

Figure 2 shows a possible ¢ = 0 contour for a system in which the 
gain of a nonelinear active device is a function of the amplitude of 
the signal. In this case, as the amplitude increases the gain decreases. 
The plane is divided into three regions for positive real values of the 
parameters a and §. The point X on the end of the M-elocus is determined 
by the maximum possible gain and the point Y by the minimum possible gain. 
If the oscillation should start with initial conditions placing the Me- 
point in the region C, the amplitude will decrease, since the oscillation 
is stable, thus causing the gain to increase and the M-point moves toward 
the ¢ = 0 contour at a. If the initial conditions are such as to put the 
Mepoint in the region B, the oscillation is unstable and will grow in 
amplitude, causing the gain to decrease and the Mepoint to move down the 
locus to the intersection with the % = 0 curve at a. Thus, regardless 
of whether the oscillation starts with M in region A or B, the amplitude 
will change until a steady state condition of a stable limit cycle is 
produced with amplitude and frequency which characterize the point of ine 
tersection at a. Should the amplitude be such that the M-point is in re- 
gion C, the region is stable, the amplitude will decrease, the gain will 
increase and the Mepoint will move up toward X and away from the inter- 
section at b. The intersection at b is also a limit cycle, but it is 
clear that no oscillation can ever be established at that frequency in 
the system described since the tendency at b is either to have the os- 
cillation collapse or move to the stable limit cycle conditions of b. 

It is seen from the above analysis, that the steady state conditions 
are such that the self-excited oscillation may be characterized by the 


intersection of the M-locus with the ¢ = 0 contour. But this, in itself 
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is equivalent to assuming that an oscillator may be analyzed by consid- 
ering its oscillation to be characterized by a pair of poles on the jw 
axis of the s-plane. Hence the application of the parameter plane meth- 
od justifies the assumption made in earlier analyses and provides, in 
terms of non-linear theory an insight into the operation of the oscil- 
lator. It should be noted that the earlier analyses were restricted to 
the steady state oscillation whereas the parameter planectheory has suf- 
ficient information to describe the transient period after turning-on 
during which the oscillation is moving to the stable limit cycle. If 
the duration of this transient or the excursions from steady state 
brought about by changes in system parameters are of interest, expres- 


sions for these changes have been developed.[ 20] 


Oscillators as Feedback Systems 
In general, any oscillator that provides continuous self-regulation 
can be considered as a feedback system. The diagram in Figure 3 shows 
such a general system in which the forward element is considered to be 
an active element in which the gain is a function of the amplitude of 
the signal. The general characteristic equation of such a system from 


block diagram algebra is 
1 - BH(s) = 0 (5) 


From this simple relationship the Barkhausen criteria for oscillation 


are immediately apparent, i.e. 


lsH(s)| = 1 (6) 


The loop gain must be unity and the total phase shift around the loop 


must be some integral multiple of 2m. 
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Figure 3 


Block Diagram of Oscillator 
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In many practical oscillators, the active element is a single stage 
vacuum tube or transistor which provides an inherent ft phase shift. In 
such a case the feedback network will provide an odd integral multiple 
of t phase shift. 

The feedback network itself may take on many forms and topologies 
and there is a fundamental problem in oscillator design to choose a suit- 
able feedback network topology and to choose parameters to suit specifi- 
cations. The parameter plane method of analysis and comparison of pro- 
perties of different topologies and a sensitivity analysis of different 
topologies provides a sound basis for choice. 

It is convenient to consider oscillators to be divided into two 
classes, those incorporating resonant circuits as the frequency deter- 
mining feedback networks and those incorporating only resistive and 
capacitive elements in the feedback circuit. The latter are important 
in the development of integrated circuits where an inductive element 
cannot be manufactured in the circuit with the ease that a resistive or 
capacitive element can. In this thesis, RC type oscillators are con- 
sidered in general with an idealized active element considered to be a 
saturating amplifier with infinite input impedance and zero output im- 
pedance, or an amplifier with deadzone. The LC oscillators considered 
are analyzed on the basis of the Norton equivalent circuit of a vacuum 
tube active device. Extension to transistorized circuits is a matter 
only of substituting a suitable model and redetermining the character- 
istic equations of the oscillators. The analysis is exactly the same 


as that to be developed. 
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4, The Sensitivity Problem in Oscillator Design. 

In practical use, oscillators are required to produce and maintain 
a waveform, often sinusoidal, of fixed frequency and magnitude. If the 
application requires that the frequency and/or amplitude be variable, 
the variation must be accurately controlled. However, there are many 
parameters in an oscillating system which are subject to variation in 
a non-deterministic manner. The reduction of sensitivity of amplitude 
and frequency to these uncontrolled changes is a basic consideration. 
The sensitivity analysis of oscillators from the parametric equations 
from the parameter plane method provides a ready comparison of topolo- 
gies and the effect of parameter values on sensitivity to parameter 
changes. In the following analysis, the effect of variation of para- 
meters is the primary consideration. It is recognized that the stabil- 
ity of the oscillator signal also depends on phase shifts in the ampli- 
fiers and changes in harmonic content of the signal. In many practical 
circuits, additional feedback compensation is made to make the ampli- 
tude relatively insensitive to frequency changes. This adaptation is 
not considered. However, it will be noted that under some circumstances 
certain unadapted systems have an inherently low amplitude sensitivity 


to frequency variation. 


Sensitivity Definition 
Sensitivity is a measure of the change in an operating character- 
istic or representation of a system with respect to the variation of a 
system parameter or parameters. The classical sensitivity problem deals 
with small parameter changes and the sensitivity operator is differen- 
tial in nature. Bodefl] is responsible for the first sensitivity defin- 


ition. Latterly, several definitions have been given and each has its 
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own utility in the context of suitable problems. 


Single Parameter Sensitivity 


The basic single parameter definition of sensitivity is 


gi change in T 


x” change in x 

The change in either T or x can be a relative change, i.e. a normalized 
or percentage change, or an absolute change. In the measure of changes 
in parameter values, when the parameters are the values of lumped net- 
work elements with percentage tolerances, to speak of a relative change 
is useful. In the definition of a sensitivity function which is the 
measure of the movement of a pole or a zero of a transfer function, the 
absolute value of the change is of significance. In the case of measur- 
ing the sensitivity of frequency of oscillators to changes in parameters 
of the oscillating system, the relative change of both frequency and 
parameter is of interest. As the problem is developed, it may be deter- 
mined that a particular sensitivity definition results in a simplifica- 
tion of expression or calculation which favors one definition or another. 


In this thesis, the Horowitz[4] definition is used as follows: 
S = aa 46 dur (8) 


This definition will be shown to lead to some simplifications in the 
calculation of frequency sensitivity of oscillators. The inherent sin- 
gularity at the null value of the parameter should be noted, but it 
causes no difficulty if it is kept in mind that this definition is a 
measure of a percentage change in T due to a percentage change in the 


parameter x. 


Multiparameter Sensitivity 


In many cases, it is highly likely that no single parameter will 
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cause all of the changes in the function of interest. A single measure 
of the cumulative effect of the variation of several parameters is re- 
quired. Goldstein and Kuof15] extended the single parameter definition 
of sensitivity into the multiparameter case. The total logarithmic dif- 
ferential can be written: 


0 alert 
ore 7) 


dinT=- 


By considering the set of fractional parameter increments as a vector, 


dY = (d ln x,, d ln x., ~, d ln x)» and setting y, = In x,, the above 


1? 2° 
total differential can be written 


dinT= V{ln T(y)> Yor ™s ne dy (10) 


If we fit this expression to our single parameter definition, the multi- 
7 = ; 
parameter sensitivity, S', 18 given by 


Si=2vlnT (11) 


The magnitude of this sensitivity vector can be expressed by the norm. 
Still this multiparameter definition leaves something to be desired 
for a general definition. In the manufacture of integrated circuits, 
for instance, all of the components are laid down at the same time under 
identical physical conditions. The value of the components can only be 
realized approximately, but, since they are all laid down together a 
high degree of correlation is expected. In this case, it might well be 
that the sensitivity with respect to some parameter is of opposite sign 
with respect to another. Indeed, good design would take maximum advan- 
tage of such a situation to reduce the variance in overall performance 
criteria for the circuit. In this case, it is not reasonable to take a 
sensitivity measure which adds the squares of the individual sensitiv- 


ities as is the case where the norm of the sensitivity vector is taken 


as the multiparameter sensitivity. Leef[26] and Biswas[3l] have recog- 
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nized this point and in Lee's concept of "Sensitivity Groups" and in 
Biswas’ inclusion of correlation terms in the definition, a more rea- 
sonable sensitivity measure for the multiparameter case is obtained. 

In the case of networks composed of individual components rather than 
an integrated circuit, the designer might safely assume independence of 
the parameter variations for which the Goldstein and Kuo multiparameter 


sensitivity definition is adequate. 


Frequency Sensitivity 

The parameter plane representation of an oscillator has been dis- 
cussed. With a definition of single parameter sensitivity in hand, it 
is now possible to see how the concept of frequency sensitivity may be 
expressed graphically in the parameter plane and algebraically from the 
parametric equations and from the characteristic equation. Figure 4 
shows a typical ¢ = 0 contour from a parameter plane drawn for an oscil- 
lator. From the previous analysis, it was seen that the intersection 
of the Melocus of the non-linearity gives the frequency of the limit 
cycle. The amplitude of the limit cycle may be given if a known rela- 
tionship exists between the amplitude and the defined non-linear para- 
meter. It is seen that every point on the % = O contour is a potential 
limit cycle point and the frequency of the limit cycle is given implice- 
itly by the parametric equation of the contour. Since the contour is 
described by the parametric equations, a@(w) and 8(w), it is possible to 


give the slope of the graph as follows: 


os = 2 ow (12) 


from which it is seen that 


04 
Boa = 8 aww ox (13) 
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M=-locus 


BIS 


2 Wy 


Figure 4 


Components of Typical © = QO Contour 


27 


or 


_ gh ow 
= = 8 Sy 


(14) 
Hence it is seen that, if B is the parameter assigned to the variable 
gain in the oscillator, due to the non-linearity, the frequency sensi- 
tivity with respect to the variable parameter, a, is given algebraically 
from the logarithmic derivative of the parametric equation @ = a(w) or 
from the slope of the © = O contour itself it there is a concomitant 
knowledge of oe 

Since the parametric equations of the { a O contour are ultimately 
determined from the characteristic equation, it will be of some utility 
to determine a relationship between frequency sensitivity with respect 


to variation of the parameters and the coefficients of the characteristic 


equation. 


Second Order System 


For the general second order system, 


2 
ays + a,s + ag = 0 (15) 


Substituting Ss = jw, it is seen that 
a 
oe send a. 0 (16) 
a 1 
a 
This gives the characteristic frequency and conditions for oscillation 


(two roots on the jw axis). The frequency sensitivity with respect to 


any parameter x, Sv is then determined from 


2 Inwae ln ay ~ in a, (17) 


Sean . 282 
5. = 5(S,,0 Sy ) (18) 
Third Order System 


For a general third order system 


28 


3 2 
as + as + a,s + ag = 0 (19) 


Substituting s = jw, it is seen that 


a a 
nie =, me (20) 


from which it follows that 


a pia e3 yn *6 a9 
Ss 2 5(8. eee ie) (21) 


It must be borne in mind that the basis of this definition in this con- 
text is the argument that there is a selfeadjusting parameter in the 
System causing the steady-state M-point to follow the @ = 0 curve. Since 
this parameter will be included in some or all of the coefficients of the 
characteristic equation, the partial derivative of a coefficient with re- 
spect to a parameter must also include the partial derivative of the 


adaptive parameter in accordance with the chain rule. 


Amplitude Sensitivity 

Since the effective gain of the non-linear active element in this 
discussion of oscillators is considered to be a monotonic function of 
amplitude of the signal, and since the effective gain is specified by 
the intersection of the M-locus with the © =» O curve, it is a temptation 
to speak generally of the amplitude sensitivity. However, the whole 
concept of amplitude with respect to oscillations in non-linear systems 
needs better definition. From the chain rule it is possible to show 


that if the effective gain is a function of amplitude only, 


B A 
s = s.-S" (22) 


A. . , Meg : ' 
where so is the amplitude sensitivity with respect to frequency. sf is 
derived from the parametric equations of the ¢ = O curve and on depends 


entirely on the relation of effective gain to amplitude. It is in the 
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relationship between amplitude and effective gain that there is diffi- 
culty. If there is sufficient low pass filtering in the feedback net- 
work, the describing function is a reasonable relation to use for the 
non-linear gain. In many cases, in practical oscillators this is not 
the case. This statement is equivalent to saying that, unless other 
steps are taken, the waveform produced cannot be nearly sinusoidal but 
will have harmonic distortion to some, possibly unacceptable degree. In 
order to reduce harmonic distortion in such cases, it is possible to ar- 
range that the linear gain portion of the non-linear gain curve is just 
high enough to guarantee that the end of the M-locus lies just inside 
the unstable region of the parameter plane. In such cases, the ampli- 
tude is limited by the signal level at which the non-linear property of 
the active device is effective to reduce the loop gain to unity. This in 
itself will cause the amplitude to be frequency insensitive but will re- 
quire a degree of sophistication in the oscillator design, if the fre- 


quency is to be variable. 


Sensitivity in Specifications 
Sensitivity specifications may be made in many forms. The most 

general specification for a fixed frequency oscillator is likely to be 
that the frequency shall not vary more than a specified percentage due 
to any reasonable cause. Due to the stochastic nature of the variation 
of parameters, this necessarily implies that the tolerance cannot be met 
by all oscillators but shall be met by some, presumably very high number. 
In addition, it is clear that the specification is essentially a multi- 
parameter specification involving the calculation of sensitivities due 

to each parameter and the probability density for each parameter. If, 


however, there is one predominant variable parameter in a system, until 
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such time as the effect of that parameter is reduced to the order of the 
variation of other parameters, itis valid to ignore the lesser effects 
and treat the problem essentially as a single parameter sensitivity 
problem. 

After a study has been made of single parameter and multiparameter 
sensitivity functions for practical oscillator circuits, it will be pos- 
sible to comment further on suitable forms for stating sensitivity speci- 


fications. 


ou 


58 Analysis of Common RC Oscillators. 

The general problem considered in this analysis is to find, on the 
basis of the comparison of parameter plane curves and sensitivity funce 
tions, the best choice of a feedback network consisting entirely of re- 
sistive and capacitive, passive devices which will yield the most stable 
frequency in the sense of minimum frequency sensitivity. The single 
parameter variation which is likely to outweigh the others is the load 
resistance. This resistance may be the input resistance of the active 
device which is likely to vary and the input resistance of the succeed- 
ing stage in any device in which the oscillator is used. In actual fact, 
the nature of the variable parameter considered is of no importance to 
the analysis since the method is identical for every parameter. However, 
conclusions drawn in the following examples are valid only for the para- 
meter values chosen and for the particular variable parameter considered. 
Each set of parameters requires its own analysis. 

The parameter plane curves are drawn for a normalized frequency 
X= RCw, thus making the analysis applicable to any frequency range. 
Arbitrarily, to provide a basis for comparison of different networks, all 
resistances are chosen equal to R other than the variable parameter, aR, 
and all capacitances equal to C. It is recognized that this may not lead 
to the best choice of parameters from other points of view. Mehta[ 16] 
has shown that this choice is not optimum to reduce sensitivity of fre- 
quency to phase change in the active element. However, this does not 
affect the method, of which this material is exemplary. 

The characteristic equation of each RC oscillator was found by the 
method of determining the transfer function of the feedback network from 
two-port network analysis and then substituting in 1 - BH(s) = 0, (5), 


the characteristic equation in general from the feedback block diagram. 
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Phase Shift Type 
The phase shift oscillator takes advantage of the inherent tt phase 
shift of a single or odd number of stages of the active element and 
provides an additional t phase shift from the feedback network to meet 
the zero phase shift criterion for oscillation. The feedback network 


may be either a high pass or a low pass configuration. 


Low Pass Configuration 

Figure 5 shows the general circuit in which the parameter of in- 
terest in the comparative analysis is marked (*). The characteristic 
equation in general is; 


3 
Ss (R) RRR, C)C,C,) 4 


2 
s (R, RAR, CC, + RRR, CC, + R,RjR4C1 6, + 


ROR4R, CC, + R,RAR,C1 0, + R,R,R,C,C,) 45 


s(R,R,C, + RoR,C, t RRC, + ROR,C, + RRC, +t 


RRC, +t RRC, + RRC, + RRC, + R,RC,) + 


(Ry . R, ar R,) te (1-8 )R, = 0 (23) 


For the particular example, setting Ry = R = Ry = R, Ry = QR, and 


9 = C, = C, and setting p = RCs, we have the reduced characteristic 


equation, 
ap? + (Sa + 1)p* + (@ + 4)pt (3+ a -a8) = 0 (24) 
By setting p= jxX, to plot the ¢ = 0 contour in the parameter plane, we 


have the parametric equations of the contour, 


4 
Cn eY" Zp 
7) P (25) 
x4 + 11X2 + 14 
a (26) 
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Figure 5 


Phase Shift Oscillator - Low Pass Configuration 





Figure 6 


Phase Shift Oscillator - High Pass Configuration 
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and calculating directly from the characteristic equation coefficients 


the frequency sensitivity with respect toq is seen to be 


S02 it 


The sign preceding 8 is a measure of the requirement that there be art 
phase shift in the active device since the network is introducing a ft 
phase shift. Figure 11 is a comparative plot of the ¢ # O contours in 
the phase plane for all the examples of RC oscillators considered. In 
order to compare the curves of the phase shift types with the others, the 


curve is plotted for \B|. 


High Pass Configuration 
This is the most common phase shift oscillator configuration. The 


network is shown in Figure 6. The general characteristic equation is: 


3 
s“[ (1-8 )R,RR,C,C,C,] + 


2 
S (R, RCC, + RoR4C50, + RRC, C, + 


R)R,C)C, * R,)RoC5C, * R,R,C,C4) + 


s(R,C + RC, + R{Ca cfr RoC) + la» 0 (28) 


i ZZ i | 2 


OK we we obtain the reduced characteristic 


equation for the particular problem as follows 


(w-0B8 )p> + (3+3x)p* + (4ta)p + 1 = 0 (29) 
Once again the third order method is used to obtain the frequency sensi- 
tivity from the coefficients of the characteristic equation. One can 
use the coefficients of the odd or even powers, but, in this case, a 
prior survey would show that the odd coefficients are functions of both 
aw and 6B, whereas the coefficients of the even terms are much simpler. 


From the general expression, 


ob) 


ear. a2 
er) Ce ie, (30) 


a Ra 
and the fact that a O, it is seen that 


l_a fo 
=~ 80? =~ Tay (31) 


Setting p= j% the parametric equations of the © = O contour are: 


2 
a= +7 (32) 
12 l 
B= jTogzZt wel (33) 


where once again the sign before B indicates the ft phase shift in the 


active element and the curve is plotted for |p]. 


Other RC Oscillator Types 
In other types of RC Oscillators, the frequency of oscillation is 
determined by the zero phase shift criterion. Each of the following 
networks: has a zero phase point where the phase changes from leading to 
lagging, or vice versa. The determination of the characteristic equation, 
parametric equations of the ( = 0 contour and sensitivity of frequency 


to a single parameter change proceed as before. 


Wien Network 
The network is shown in Figure 7. The general characteristic equa- 


tion is: 
s“(R)R,C,C,) + s[R,C, + RC, + (1-8)RQC CJ o+ 12 0 (34) 


Setting Ri=R; Rar, C#C.ac we obtain the reduced characteristic equa- 
tion 
2s ll Ocal sues Sec agi Ey (35) 


from which, setting p= jX, we find 
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Figure 7 


Wien Network 


Figure 8 
Bridged "T'' Network = Type 1 


Ey 


and gf = ° ; 


(36) 


(37) 


(38) 


It is worth noting that the network is extremely simple but at the 


cost of having the sensitivity of frequency with respect to any parameter 
change being a constant value higher than that in either of the two phase 


shift configurations thus far considered. For the sake of comparison, 


the frequency sensitivities are plotted in Figure 12. 


Bridged "T'' Network - Type l. 


Since there is a symmetry between resistive and capacitive elements 


in Figures 8 and 9, the configuration of Figure 8 is arbitrarily called 


type 1 for the sake of identification. The general characteristic equa- 


tion is: 
s7[ (1-8 )R,R,R.C,C,] + 
ji i ie ga? 


s[ (1-8 )(R,R4C + R.R.C_) + RRC, + R,R,C, + R,RAC,] + 


1 23 2 32 | beta P22 
R, + (1-8)R, = 0 (39) 
Substituting as before, the reduced characteristic equation is 
(w-08 )p” + (2+3a-208)p + (1te-08) = 0 (40) 


from which it is determined that the parametric equations of the @ = 0 


contour are 


2 


(41) 
14207 


and the frequency sensitivity with respect toad is: 
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Figure 9 
Bridged "T'' Network = Type 2 


Figure 10 


Twin "T'' Network 
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2+01 (43) 


Bridged "T" Network =- Type 2. 
This network is shown in Figure 9. The general characteristic equa- 


tion is: 


sf (1-6 )R1R,R4C,C.] + 


s[ (1-8) (R,R,C, Zs RRC) 7 R, RAC, + R,R,C,] + 


R, + Ro a (1-6 )R, = 0 (44) 
Substituting as before, the reduced characteristic equation is: 
(v-o8 )p* + (1+37-208)p + (2ta-~aB) = 0 (45) 


from which it is determined that the parametric equations of the & = 0 


contour are 


2 
ae st es (46) 
2 
= =a» (48) 
ane of = (atl) @-3) 


Twin "T'' Network 
This network is shown in Figure 10. The general characteristic 


equation is: 


s7T (1-8 )R, RR4R,C) CoC] + 


2 
S [ (1B) (R)R,R,C,C, + R,R4R,C,C.) + R, RoR, CoC, ct: Ri RR301 6, + 


R)R,R4C,C, + R)R4R,C,C, + R,R4R,C C3] “ 


s{ (1-8) (RRC) + R4R,C.) + R,R4C. + R,R,C, + RRC, + R, RC) + 


R,R,C, + R,R,C, + R)R,C, + R,R,C,] + 


40 
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CR, + R, (1-8 )R, J = 0 (49) 
Substituting as before, the reduced characteristic equation is: 


3 2 
(a-a8 )po + (2+5a-208 )p + (St5a-208 )p + (2tv-08) = 0 (50) 
from which it is determined that the parametric equations of the f= 0 


contour are 


~ 2. 
am we (51) 
64 9X4-2172- 


and “i = = (so) (53) 


6 
where S12 = - = x Boge a (54) 


It is worth noting that this last analysis maybe about as difficult 
an analysis as one might consider making using algebraic operations with- 
out the aid of a computer. If the computer is available, the generalized 


sensitivity equations of Siljak and Stoic[21] may be used. 


Stability Analysis of Limit Cycles 

As noted in the theory of self-excited non-linear oscillations in 
the parameter plane, the stable regions of the parameter plane as parti- 
tioned by the % = O contour, must be examined. If we assume that the non- 
linearity is such that the parameter, B, is a monotone decreasing func- 
tion of the amplitude of the signal, it is quite easy to take test points 
above and below the & = 0 curve to determine the stability of the limit 
cycles. On this it can be shown that both phase shift networks and the 
Wien Network give rise to stable limit cycles and both bridged "T" types 
and the twin "T"’ network give rise to unstable limit cycles. If a non- 
linearity is chosen in which the gain increases as the amplitude increases, 


such as a dead-zone type, the latter three networks give rise to stable 
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limit cycles. Appendix A describes a basic problem in the parameter 


plane method arising from this analysis 


Example 

Consider the details of a possible design analysis in the parameter 
plane based on the curves developed thus far. An oscillator is to be 
built using an amplifier with a linear portion gain of 70. A phase 
shift, high pass configuration has been chosen. It is desired to mini- 
mize sensitivity to the network terminal resistance. (See Figure 6.) 
What should @ be and what is the effect on the sensitivities with respect 
to the other parameters if aw is changed from an arbitrary reference at 
Y= 1? 

The sensitivities are calculated from the coefficients of the char- 


acteristic equation. It is seen that 


sos - 3902 (55) 


gO _ oO (56) 


The sensitivities at aw = l1 are: 


cs . _ - OLS se . - O63395 s — - 0.250 
R, R, R, 
(57) 
W W W 
SC, = ~ 02250 Sc, - - 0.333 5c, ee 1 IG 


The maximum possible sensitivity with respect to any single parameter is 
0.5 (see Figure 12). From Figure 12, sensitivity with respect to R, is 

minimized by making @ as small as possible. However, references; to Fig- 
11 shows »dmmediately, that B increases as a is decreased, hence the 


linear portion gain of 70 is an upper bound on 8, implying a lower bound 


ond. A suitable value for a is 0.3 for which B = 55. The sensitivities 
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ata = 0.3 are: 


Sp RAG SR eal S. tes Opell 
1 2 3 
(58) 
Sie ei0295 SY eres so Sen 0Ns72 
C, Gy C. 


It should be noted that although the required minimization has been ef- 
fected within the constraints of the problem, the sensitivity with ree- 
spect to three others has been increased! This, apparently, is yet 
another application of the universal principle that you can't get some=- 
thing for nothing: 

Reference to Figures 11 and 12 also shows immediately that, in the 
problem above, second thought should be given to the choice of a high 
pass phase shift network. The low pass phase shift network would have 
provided minimization of sensitivity to R, much lower than the 0.115 
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figure reached for any value of aw greater than 2.5. The value of R, is 
very likely to be determined by the amplifier specified. If it is suit- 
able to choose R less than Ra /a for the required frequency, then a change 
to the other phase shift configuration is recommended. If this choice 

is not suitable (the required C may be very large) then it may well be 


that the phase shift high pass configuration may be the only one capable 


of handling the low input impedance with a low sensitivity. 


General Observations 
A scrutiny of Figures 11 and 12 reveals many of the characteristics 
of the RC Oscillators studied. For those oscillators with a portion of 
the © = O contour almost parallel to the q@ axis, it is seen that dB/d& 
is very nearly zero. Since the frequency sensitivity with respect to 
B 


w is not zero in these portions, it must be that SO is very nearly zero. 


This implies that the amplitude sensitivity is very nearly zero in these 


ae, 


regions, regardless of the relationship between § and amplitude. Hence, 
to minimize amplitude sensitivity to change of frequency in a variable 
frequency oscillator, it is worthwhile to design in a flat portion of 
the © = O contour. 

It is immediately seen from the parameter plane and sensitivity 
curves which configurations can meet an a priori single parameter sensi- 
tivity, and the required amplifier gain is also shown. The Melocus for 
the non-linear amplifier indicates whether it is feasible to consider a 
configuration, since there may not be an intersection with the f = 0 
contour, or if there is such an intersection, the indicated limit cycle 
may not be stable, or the sensitivity specified may require an impossibly 
high gain. A set of curves can be drawn for any parameter of interest 


using suitable preselected values of the other parameters. 


Multiparameter Sensitivity 

The comparisons in Figure 12 were based on the single parameter 
sensitivity of frequency with reSpect to a terminal resistance in the 
feedback network. Just so long as the variation of this terminal resis- 
tance is the predominant variation of parameters in the oscillator, then 
it is valid to make comparisons of configurations based on this single 
sensitivity. However, it is apparent that, when there is a likelihood 
of several parameters varying with an equal effect, then the multipara- 
meter sensitivity measure must be taken. Provided the parameter varia- 
tions are uncorrelated, the Goldstein and Kuo definition is satisfactory. 
Since this definition of sensitivity is a vector sensitivity, the maxi- 
mum sensitivity is given by the norm of the sensitivity vector. From the 
definition 


d ln T # V(ln T)°(d ln x,, », d ln x) (10) 


Ge 
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where = VOR MD Sing, “17 Oink “a ee 
L mer Te .2 4 
and SS | ase Ji) 2 (59) 


On this basis, it is possible to compare the multiparameter sensitivi- 
ties of the two types of phase shift networks considered in the discus- 
sion. 

It is assumed that the system parameters have been chosen in ace 
cordance with some other performance criteria and it is desired to com- 
pare the multiparameter sensitivities of the two phase shift types and 
Wien bridge type of feedback network. Arbitrarily, for the sake of exe 
ample, it is assumed that all resistances are equal to R and all capace 
itances are equal to C. This does not fix the frequency, hence the 
problem is not completely abstracted from reality: The reduced charac- 
teristic equation is determined for each parameter considered in turn to 
be variable or, to be even more general, the general single parameter 
sensitivities are determined from the general characteristic equation 
for all system parameters. The norm of the sensitivity vector is de=- 
termined from the components (single parameter sensitivities). Table l 
shows the coefficients of the reduced characteristic equations for each 
of the three types considered. It should be noted that this is not a 
minimum sensitivity but a demonstration of the feasibility of getting a 
multiparameter Sensitivity index and the fact that it can be used as a 
basis of comparison of feedback networks. 

From Table 1, it can be deduced that for a = 1, the multiparameter 
sensitivity of the phase shift, low pass configuration network is less 
than either of those of the other two types considered. This is similar 


to the trend already apparent in the single parameter sensitivity to the 
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Table l 


Determination of Multiparameter Sensitivity 


Phase Shift Oscillator-Low Pass Configuration 


Parameter 
Ry 
2 
3 
Ry 
On 
ey 
3 

Phase Shift 
ae: 
R, 
R, a 
C, or 
C. o 
C. 4 


a 


eg 


ha + 2 


4a + 2 


Ga + 2 


om © & 


4a 


3 


Oscillator-High Pass 


- a6 


- ap 


- ap 


- ap 


cope eae | 


4a + 2 


Sopwnioon3 


Wien Network Oscillator 


For all parameters, S.. 


201 


20 


2a 


ail *0 
+a tee - 6 
+6 a+3e-8 
+6 @+to-p 
+4 a+ 3 = ap 
+ 7 4-8 

+ 6 4-68 
= 4-68 
Configuration 

+ 3 1 

+ 3 1 

+ 4 1 

+ 4 1 

‘le 1 

+3 1! 
- - Hence, S = 1.000 
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i = I — i — 
Diy Slt ple Ble sit Blu 


S= 0.833 


terminal resistance for @ = 1 shown in Figure 12. Once again, curves 
might be drawn for comparison purposes, but it is recommended that a gen- 
eral purpose computer be used if one is available since the calculation 
of one point on the curve represents considerable time spent in calcula= 
tions. 

It should be noted that several assumptions have been made in de= 
termining this sensitivity index for many varying parameters. It has 
been assumed that each of the parameters is independent, 4 questionable 
assumption if one is considering the effect of ambient temperature on 
frequency stability since it is highly likely that all similar components 
would suffer similar variation when subjected to the same variation in 
temperature. It has been assumed that the variation (relative or incre=- 
mental) for each parameter is the same. In the case of the stochastic 
nature of the variation of parameters, this assumption is equivalent to 
assuming that all parameters have the same variance, valid, perhaps for 
the passive elements in the network but questionable when one considers 
the terminal resistance of the network. If each of the incremental vare 
iations of parameters is known, a better assessment of the incremental 
change in frequency to be expected can be made. However, to keep the 
work of determining system sensitivity to a minimum it should be noted 
again that the predominance of any single parameter, renders variations 
due to other parameters secondary and the problem might be considered a 


single parameter sensitivity problem with much time saved. 


Algebraic Determination of Parameters 
Thus far, the discussion has been based on the assumption that the 
oscillator design has been substantially completed for several configure 


ations and then the sensitivity analysis is applied to determine which 
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design is best suited to the requirements. In each case it is possible 
to determine which configuration represents the minimum sensitivity to a 
parameter of interest but there is no indication that this sensitivity 
is the minimum for that system. Indeed, it might be possible to find a 
minimum sensitivity set of parameters for one of the configurations re- 
jected that might make it more suitable. Lest it be considered that 
this is a simple min-max problem, note.that there are physical con- 
Straints on the solution of the minimum sensitivity problem, such as 

the fact that all parameters must be real, positive numbers; capacitors 
larger than about 100 mfd and resistors larger than 22 megohms are 
specialty items; amplifier gains for simple one or two stage active de- 
vices are seldom over 100; the nature of the sensitivity functions al- 
ready analyzed show functions which do not have minima for finite values 
of the argument. In some very simple cases, it is possible to derive 
general sensitivity functions. By finding the conditions for a minimun, 
one equation of n possible equation is determined to establish the values 
of the n parameters. The frequency required is not considered as one of 
the equations in this context for RC oscillators since our entire devel- 
opment has been based on a normalized frequency which is determined by 
the choice of R and C after the design is optimized. 

One need only look at the nature of the coefficients of the general 
characteristic equations in the previous examples to realize that the 
sensitivity functions are of such a complexity that it would take several 
hours to set up the min-max problem. It appears that an algebraic deter- 
mination to the choice of parameters for a minimum of sensitivity is a 


¢ 


problem well suited to a computer solution. 
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Figure 13 


Tuned-Plate Oscillator 
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Figure 14 


Tuned-Plate Oscillator = Equivalent Circuit 
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6. Analysis of Common L-C Oscillators. 

The elements of the application of sensitivity analysis to oscil- 
lator design have already been developed and demonstrated for several 
RC oscillators. However, for the sake of completeness and to point out 
some of the essential differences the class of oscillators incorporating 
resonant circuits is considered. The active device in these oscillators 
is the non-linear device once again but in this case there is not the 
same idealization. The active device is considered to be a vacuum tube 
and the Norton equivalent (current source) is used in the analysis. The 
justification for using the linear analysis on an essentially non-linear 
problem is as before; the oscillation is a stable limit cycle and can be 
characterized in the parameter plane by the steady-state position of the 
M-point moving on the {= O contour. The characteristic equations aris- 
ing are second and third order allowing the application of the simplified 


calculation of sensitivities from the characteristic equations. 


Tuned Plate Oscillator 

A resonant tank circuit in the plate circuit is inductively coupled 
to the grid circuit to provide the necessary positive feedback to sus- 
tain oscillation. The circuit diagram and equivalent network are shown 


in Figures 13 and 14. The general characteristic equation is: 


2 
Ss (LC) + s(Liy, + Cr+ gM) + (1 + ry) = 0 (60) 


where y= (61) 


Ona 


a Me Lo 


gl 

p *&L 

It is worth noting that the characteristic equation coefficients do not 

9 but only on the mytual inductance M. In order for the cire- 
D 


cuit to oscillate, the mutual inductance, M, must be negative, i.e. there 


depend on L 


is a right and a wrong way to hook up the coupling transformer. This is 
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Figure 15 


TunedeGrid Oscillator 
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Figure 16 


Tuned-Grid Oscillator = Equivalent Circuit 
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seen from the coefficient of the first order term which must vanish for 
oscillation and cannot do so unless the mutual inductance is negative. 
This might also be stated that the vacuum tube introduces af phase 
shift which must be compensated by aft shift through the coupling trans- 
former in order to satisfy the unity loop gain, zero phase Barkhausen 


criteria for oscillation. 


Tuned Grid Oscillator 

A resonant tank circuit in the grid circuit inductively coupled to 
the coil in the plate circuit provides the required feedback. The cir- 
cuit diagram and equivalent circuit are shown in Figures 15 and 16. The 


general characteristic equation is 


3 Zz Z 
s"(C(L,L, - My] + s*[c(L,ry, + L,)] + 
Ser Ly, + gM) + l= 0 (62) 


Colpitts Oscillator 

This oscillator and the Hartley oscillator are very similar in that 
the plate signal is coupled back into the grid by a balanced resonant 
circuit. The Hartley is more difficult to analyze due to the inevitable 
resistance of the inductors and also due to the mutual inductance be- 
tween the inductors, which does not exist in the Colpitts case. The 
Colpitts circuit and equivalent circuit are shown in Figures 17 and 18. 


The general characteristic equation is: 


3 2 
s-(LC,C,) + s (LC,y, + rC,C,) + 
s(rC,y, + Cc, = Co) + co + g = 0 (63) 


Hartley Oscillator 


The circuit and equivalent diagrams are shown in Figures 19 and 
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Figure 1/7 


Colpitts Oscillator 


Figure 18 


Colpitts Oscillator ~- Equivalent Circuit 
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20. The general characteristic equation is: 


3 2 2 
S (Ce LL. + cy LL, Ce M Cy M ) + 


2 
S (CL) rog. + CL rié,+ cLyroy,, + CLOT VY, + CL, + CL. + 2CM) + 


2 ] 2 


s(Cg rir, + Cy oT 16 + Loy, + Cr, + Cr, = gM) + (roy,t 1) = oO 
(64) 
Sensitivity Analysis 
At first sight, the sensitivity analysis appears rather straight 
forward. The equations are second or third order and the sensitivities 
with respect to any parameter can be determined from the coefficients. 
However, the use of the Norton equivalent of the vacuum tube is a step 
closer to reality in this analysis and there are now two parameters 
which are not independent. To make matters a little worse, these two 
parameters are also the non-linear gn and ry of the vacuum tube. A crite 
ical analysis of the vacuum tube characteristics will show which of the 
two (g., or rt) is the least linear and presumably the most effective in 
stabilizing the amplitude of oscillation. This parameter should be 
chosen as the nonelinearity in the parameter plane representation. In 
the case where the second nonelinear parameter is a parameter of inter- 
est, there may be a method of constructing the Melocus to find the ine 
tersection with the © = 0 contour, but only if there is a known rela- 
tion between the parameters and the signal amplitude. Little error is 
introduced by assuming the least effective nonelinear parameter to be 
constant for the sake of the analysis. 
Another important difference presents itself in that the RC nete 
works lent themselves to a very convenient normalization which made a 


single analysis applicable to oscillators at any frequency. In the case 
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Figure 19 
Hartley Oscillator 
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Figure 20 


Hartley Oscillator + Equivalent Circuit 
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of the LC oscillators, the characteristic equations do not lend them- 
selves to such a simplification. Even so, sensitivity analysis can 


lead to some immediate conclusions concerning these oscillators. 


Example 
The element most likely to vary in these oscillators is the induc- 
tance. If the tuned plate and tuned grid oscillators are compared, it 


is seen that for the tuned plate oscillator 


S =-3 (65) 


whereas, for the tuned grid oscillator, 


Ww ih LiCryo 
I Pam a emer (66) 
Ly 2 C(L)ry, + L,) 
Ww 1 Lac 
| ee ed (67) 
L, 2 C(Lyry, + L,) 
If the variations of Li and L, in the tuned grid oscillator are uncor- 
related, the multiparameter sensitivity with respect to Ly and L, will 
be less than + . However, if they are correlated the multiparameter sen- 
sitivity will be equal oe . This is derived from the Goldstein and 


2 


Kuo definition of sensitivity. Hence, it appears, that there is likely 
to be a smaller frequency sensitivity due to changes in inductance in 


the case of the tuned grid circuit. 
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de Conclusion. 

Parameter plane analysis and sensitivity analysis for the frequency 
characteristics of oscillators can be carried out quite readily for syse= 
tems with second and third order characteristic equations, without re= 
sort to lengthy and generalized parameter plane expressions. Comparison 
of the © = O contours in the parameter plane gives a useful basis of com- 
parison of alternative structures and reveals the frequency of oscilla- 
tion and indicates areas of amplitude sensitivity. Unless there is a 
known relationship between amplitude and the nonelinear gain or amplifi- 
cation factors, the amplitude cannot be predicted. Only in cases where 
there is sufficient low pass filtering is the describing function an ade- 
quate relationship between gain and amplitude. In most oscillator nete- 
works, this filtering requirement is not met. 

The analysis is easiest to carry out if there is one predominately 
varying parameter, i.e. it is a single parameter sensitivity problem. A 
multiparameter sensitivity index defined by Goldstein and Kuo is an ef- 
fective comparison but requires assumptions concerning independence and 
similar variances for all parameters, to be meaningful. 

Each oscillator design problem has its own constraints. In some 
cases, minimum sensitivity to parameter variation is desired, in another, 
for a variable frequency oscillator, a wide range of frequencies with 
linear parameter variation is required. Consequently, no sweeping gener= 
alities can be made as to which oscillator circuit is best from either 
point of view. However, within the context of each design problem, there 
is an optimum configuration for frequency sensitivity to parameter variae 
tions and such a configuration may be determined by the methods developed, 
provided that the system characteristic equations are not of higher order 


than the third. Most of the well known oscillator circuits meet this 
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criterion. The characteristic equations of ten common oscillator cir- 
cuits are developed and the method is easily adapted to variations of 
these circuits. 

The method employed suggests applications to the frequency charac- 
teristics of active and passive filters in that the poles on the jw axis 
(or the points on the { = 0 contour) might equally well be the poles in- 
troduced by a sinusoidal driving function. 

The worked example shows that the single parameter sensitivities are 
functions of the system parameters. If the sensitivity with respect to 
one parameter is reduced, it is highly likely that the sensitivity with 
respect to another has been increased. This possibility must be consid- 


ered in any application of the methods of this paper. 
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APPENDIX A 


A BASIC PROBLEM ARISING FROM THE 
STABILITY ANALYSIS OF LIMIT CYCLES 


IN THE PARAMETER PLANE 


Introduction 

In the complex frequency plane, a fundamental criterion of absolute 
stability is that there shall be no roots of the System characteristic 
equation lying in the right half plane. The imaginary axis, or 0 = 0, 
line in the complex plane is the boundary between the regions of stabil- 
ity and instability. This boundary can be mapped into the parameter 
plane as a (@ = 0 contour. In parameter plane theory, this contour is 
considered to be the boundary between regions of stable solutions and 
unstable regions. However, an example arising out of the analysis of 
the stability of limit cycles in the oscillator systems, in this thesis, 
indicates that this is not the case and that there are other boundaries, 
Other than the 0 = O boundary in the parameter plane to mark off regions 


of stability and instability. 


The Particular Problem 
The reduced characteristic equation for an RC oscillator using a 


bridged "T" network as a feedback network is given by 


(u-0B8 )p* + (2430-208 )p + (1ta-08) = 0 (40) 
The first step is to find the parametric equations of the 0 = O contour 
in the a, $8 parameter plane. By the method employed in the thesis, the 
argument is offered that f = 0 for all cases when the coefficient of the 
first order term in (40) vanishes. Hence the non-parametric equation of 


the contour is given by: 
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from which it is seen that 


1,3 
pects (69) 


This is the equation of a rectangular hyperbola with asymptotes @ = 0 
and § = 1.5. Substituting s = jX im the characteristic equation and 
setting the real and imaginary parts equal to zero has the same conse- 
quence. This is by way of insisting that the solutions be a pair of 


purely imaginary complex roots. This leads to the equations: 


(a-08 )X* = (l+0-08 ) (70) 
(2+ 37-208) = 0 (68) 


The parametric equations of the % = 0 contour are then 


(41) 


QR 
if 
mS 
NO 


8 = dae (42) 

Lest there be some argument that this does not come from the general 
parameter plane method, the equations will be developed from the gener- 
alized parameter plane method. Table 2 shows the development of the 
necessary intermediate functions to find the parametric equations of the 
C= O contour. It is seen that the parameter plane method in general 
maps the imaginary axis into the rectangular hyperbola developed before 
from basic concepts (and a simpler method) and into a point at infinity, 
(0, -). 

If the € = O contour is the boundary between regions of stability 
and instability, a set of test points should reveal the existence of 
stable and unstable solutions in accordance with this. This is done for 


5 points, a = 10, 8B = O, 0.5, 1, 1.5, 2. The Routh=Hurwitz stability 


ol 


criterion is applied for each set of coefficients determined. This de- 


velopment is shown in Table 3. 


Conclusion 
The { = 0 curve does not provide the only boundary in the parameter 


plane between regions of stable and unstable solutions to the system 


differential equation. 
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Table 2 


Parameter Plane Analysis - Example 


a, = bat cB + hop + d, 


k ek k 
X) = Ase? XW Uy.) Xo = wa) x WU, (C) 





k b c h a 
2 1 0 -1 0 
ii 3 0 -2 2 
0 1 0 -1 1 
B ie 1 Baw = 4 
Cc, = 9 Cc, = 0 
D=-l Do = - 2 
See 1 fe 
i= W 5 = ow 
3 
a=zw "7 W 
b= O 
c= 0 
ieee ieee 
est f- 2w> 
3 
JI= aw -w) # 0 
-e + Je2 - 4ac= Oor-= 
a = a 
eZ 2a 
2 
2w 
wy o> Tue 
8 ae BIO] 2 + D] ~~ e BO 1.2 + D2 
ibd, a ee al Cc) Hy ot C, 
1+ 2w2 
P,=° Po aa, 
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Table 3 


Routh Stability Criterion - Example 


Substitution of points in parameter plane and application of 
Routh-Hurwitz stability criterion. For second order, the Routh array 


in general is 


ee 
ay 0 
ao 0 
0 


i.e. if there is no sign change in the sequence of coefficients the 
system is stable. Number of sign changes is number of roots in right 


half plane. @ = 10 throughout. 


B ao a) ay Roots in RHP 
0 10 32 11 0 stable 
0.5 5 22 6 0 stable 
1.0 0 12 1 transition 
1.0 = 0+ Jas a i+ 0 stable 
1.0 + O- 12 = le 1 unstable 
5 © 5 2 - 4 2 unstable 
2.0 -10 - 8 - 9 0 stable 


Figure 21 shows the @ = 0-contour and the points considered. 
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Scape 





Figure 21 
C = O Contour for Bridged "T'' Network Oscillator 
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Suggestions 

The mapping of the complex frequency plane into the parameter plane 
is such that some points map into lines. It is possible for roots to 
appear in the right half plane after leaving the left half plane at in- 
finity, hence, the point at infinity must be carefully mapped in the 
parameter plane. At the origin of the complex frequency plane it is im- 
possible to speak of the value of ¢ since ¢ is indeterminate at the ori- 
gin. This raises the possibility that the origin must be separately 
mapped into the parameter plane, to dilineate areas of different numbers 


of roots in the left half plane. 


Mapping the point at infinity yields 
(a -a8 ) = QO 
which is the pair of lines a = 0 and Bs # 1. Both of these lines are 


significant boundaries when considering the Routh criterion. 


Mapping the origin yields 
(l +a -ap) = 0 
This maps into another rectangular hyperbola asymptotic toa = O and 
B = 1. This line also yields a boundary between regions of varying num- 
bers of roots in the left half plane. 

If the Jacobian is considered in accordance with 2], it is seen 
that there is a singularity atq@ = 0. This is confirmation of the above 
boundaries. 

It is evident that one must be very careful not to generalize about 
the absolute nature of the = 0 contour as the boundary between regions 
of stability and instability. This represents a boundary only for the 
case of a pair of complex roots moving into the right half plane of the 


complex frequency plane. 
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